x=[0.00025;0.00025;0.00075;0.00075];
y=[0;2000;2000;0];
plot(x,y);
xlabel('L_{tbc} (m)');
ylabel('PDF(L_{tbc})');